Energy spectra of developed superfluid turbulence 
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Turbulence spectra in superfluids are modified by the nonlinear energy dissipation caused by the 
mutual friction between quantized vortices and the normal component of the liquid. We have found 
a new state of fully developed turbulence which occurs in some range of two Reynolds parameters 
characterizing the superfluid flow. This state displays both the Kolmogorov-Obukhov |-scaling law 
Ek oc k~^^^ and a new "3-scaling law" Ek oc k~^, each in a well-separated range of k. 

PACS numbers: 43.37.+q,47.32.Cc, 67.40.Vs, 67.57.Fg 



Superfluid consists of mutually penetrating compo- 
nents - viscous normal component and one or several 
frictionless superfluid components. This explains why 
different types of turbulent motion are possible depend- 
ing on whether the normal and the superfluid compo- 
nents move together or separately. Here, we are inter- 
ested in the most simple case when the dynamics of the 
normal component can be neglected. This occurs, for ex- 
ample, in the superfluid phases of "^He where the normal 
component is so viscous that it is practically clamped to 
the container walls. The role of the normal component 
in this case is to provide the preferred heat-bath refer- 
ence frame, where the normal component, and thus the 
heat bath, are at rest. Dissipation takes place when the 
vortices move with respect to this reference frame. Tur- 
bulence in such a superfluid component with the normal 
component at rest will be called here superfluid turbu- 
lence. 

Recent experiments in '^He-B demonstrated that the 
fate of a few vortices injected into a rapidly moving super- 
fluid depends on a dimensionless intrinsic temperature- 
dependent parameter q rather than on the flow velocity. 
At q 1, a rather sharp transition is observed between 
laminar evolution of the injected vortices and a turbulent 
many-vortex state of the whole superfluid. This adds a 
new twist to the general theory of turbulence in super- 
fluids developed by Vinen and others. Attempts to 
modify the theory in order to incorporate the new phe- 
nomenon, have been made in Refs. Q, and 0|. 

In this Letter we describe how the celebrated 
Kolmogorov-Obukhov |-law for the turbulent energy 
spectrum in normal fluid, Ek oc k^^^^, gets modifled in 
the superfluid turbulence, giving rise the much steeper 
decrease, Ei. cx k~^. 



As a starting point we utilize a coarse-grained hydro- 
dynamic equation for the superfluid dynamics with dis- 
tributed vortices. In this equation the parameter q char- 
acterizes the friction force between the superfluid and 
the normal components of the liquid, which is mediated 
by quantized vortices. According to this equation, tur- 
bulence develops only if the friction is relatively small 
compared to the inertial term, i.e. when q < 1. Here, we 
will study the case of developed turbulence which must 
occur at g ^ 1. 

An important feature of superfluid turbulence is that 
the vorticity of the superfluid component is quantized 
in terms of the elementary circulation quantum k (in 
■^He-B, K = nh/m where m is the mass of '^He atom). 
Thus, superfluid turbulence is a chaotic motion of well- 
determined and well-separated vortex fllaments 3]. Us- 
ing this as starting point we can simulate the main ingre- 
dients of classical turbulence - the chaotic dynamics of 
the vortex degrees of freedom of the liquid. However, to 
make the analogy useful for classical turbulence one must 
choose the regime described by equations of the hydrody- 
namic type valid at length-scales above the inter-vortex 
distance, the latter being a microscopic cut-off similar to 
the inter-atomic distance in conventional hydrodynam- 
ics. The coarse-grained hydrodynamic equation for the 
superfluid component is obtained from the Euler equa- 
tion for the superfluid velocity v = Vg after averaging 
over the vortex lines (see review 0): 



dv 
'dt 



{v ■ V)v + = Z3 , 



(1) 



where is the chemical potential and D describes the 
mutual friction: 

D = — q;'(v — Vn) xu; + acl)x[u;x(v — Vn)] . (2) 
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Here a; = V x v is the superfluid vorticity; cl) = w/w; Vn 
is the velocity of the normal component (which is fixed); 
a' and a are dimensionless parameters describing the mu- 
tual friction between superfluid and normal components 
of the liquid mediated by quantized vortices which trans- 
fer momenta from the superfluid to the normal subsytem. 
For the flow with vortices locally aligned with each other 
these parameters enter the reactive and dissipative forces 
acting on a vortex line as it moves with respect to the 
normal component. For vortices in fermionic systems (su- 
perfluid "^He and superconductors) such forces acting on a 
vortex were calculated by Kopnin , and they were mea- 
sured in ^He-B over a broad temperature range by Bevan 
et al. . Here we consider a' and a as phenomenological 
parameters, assuming the general case where quantized 
vortices are not aligned locally and thus the bare param- 
eters are renormalized. 

Further we shall work in the reference frame where 
Vji = 0. In this frame, the nondissipative flrst term in 
Eq. renormalizes the inertial term v x a; in the left 
hand side (LHS) of Eq. ^ by the factor 1-a'. The role of 
the Reynolds number in this hydrodynamics, i.e. relative 
magnitude of the two non-linear terms, the inertial and 
friction ones, is played by the velocity independent ratio 
of dimensionless parameters, 7?e = (1 — a') /a, which was 
denoted in Ref. .1] as Tie = 1/q. The role of the parame- 
ter q as the inverse Reynolds number was demonstrated 
in experiments of Ref. 1] , where it was shown that the 
turbulence develops only below some critical value of q 
(i.e. at q < qc ^ 1). 

Here we are interested in the region of large Reynolds 
numbers, 7?e ^ 1 (9 <C 1), where the inertial term is 
strongly dominating. In this region one expects well 
developed turbulence characterised by a Richardson- 
Kolmogorov-type cascade which is modified due to the 
non-linear dissipation. In ■^He-B, the range q ^ I occurs 
at low temperatures, where a' <^ a, q ^ a. Then the 
mutual friction term in Eq. ^ can be written as 



D = qcj X [lj X v] 



(3) 



and we finally arrived at the hydrodynamic equation (^) 
whose LHS is the usual Euler equation with the nonlinear 
term which is responsible for the energy cascade in the 
developed hydrodynamic turbulence, while the RHS con- 
tains the nonlinear dissipation term given be Eq. Q . We 
will describe this cascade in the simplest possible manner, 
using the differential form [i(tI of the energy transfer 
term in the energy budget equation (the more compli- 
cated version with second derivative was used for 
superfluid turbulence by Vinen Q): 



dEk 
dt 



dk 



(4) 



Here Ek is the one-dimensional density of the turbulent 
kinetic energy in the fc-space, defined such the the total 



energy density (in the physical space) E is given by 

E^\{\v\^)^ j dkE^. (5) 



Dissipation of energy on scale k is described by the Vk 
term in the right hand side (RHS) of Eq. Q which will 
be clarified later. The idea of [lOj is to relate Ek and 
Ek in Eq. (@J in the spirit of Kolmogorov 1941 (K41) 
dimensional reasoning: 



Ek 



2/3.-5/3 



(6) 



Here C ~ 1 is the Kolmogorov dimensionless constant. In 
the absence of dissipation, Eq. Q| immediately produces 
the stationary solution Sk — s with constant energy fiux 
e in the inertial interval of scales. Then Eq. © turns 
into the Kolmogorov-Obukhov |-law for Ek'. 



Ek 



Ce2/3fc-5/3 . 



(7) 



The goal of this Letter is to describe possible modifica- 
tion of the scaling exponents in Ek due to different, than 
in the Navier-Stokes equation, form (^J of the dissipation 
term. The energy balance equation in the differential ap- 
proximation (0J is just an adequate tool for this study: 
it is as simple as possible, but not more. More accu- 
rate integral representation for the energy transfer term 
gives exactly the same results for the scaling exponents, 
because they arise from the power counting. Clearly, ap- 
proximation ^ docs not accurately control numerical 
prefactors and the exact functional form of the possi- 
ble crossover region, but these are not important for the 
questions we aim to study here. 

Within the same level of accuracy, we can simplify the 
vectorial structure of the dissipation term D and average 
Eq. over the directions of the vorticity lj (at fixed 
direction of v) 



--q\uj\v 



-q\uj\v 



(8) 



Again, we are not bothered by the numbers and there- 
fore skipped for simplicity factor 2/3 in the last of Eqs. 

. Notice that vorticity in hydrodynamic turbulence is 
usually dominated by the fc-eddies (i.e. motions of scale 
^ 1/k) with the largest characteristic wave- vector fcmax- 
These eddies have the smallest turnover time r^-m that 
is of the order of their decorrelation time. One can show 
that the main contribution to the velocity (not vortic- 
ity) in the equation for the dissipation of the A;-eddies, 
Pfc, with intermediate wave- vectors fc, 1/R < k <ti fcmax, 
is dominated by the fc'-eddies with k' ^ k. Because 
the turnover time of these eddies Tk' ^ TVnin, we can 
think of in Eq. on time intervals of interest 
(Tmin "C T ^ Tfc/ ) as sclf-averagiug quantity because it 
is almost uncorrelated with the velocity v which can be 
treated as dynamical variable. In this study, this allows 
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us to neglect in Eq. ^ the fluctuating part of \uj\ and 
to replace |a;| by its mean value. In this approximation 
Eq. ||HJ) takes very simple form: 

Tv, T =qiJo, ujo = {H) ■ (9) 

From Eq. one easily finds that 

Vk = TEk, (10) 

and the balance equation Q in the steady state finally 
takes the form: 



_ „p 2/3,-5/3 



(11) 



in which for the simplicity we put C = 1, because in our 
simple approach we are not controlling numbers of the 
order of unity. 

Let us analyze the solutions of Eq. that arise in 
presence of a fixed energy influx = e_|_ at A: = 1/ R into 
the turbulent system. Hereafter R is the outer scale of 
turbulence that of the order of the radius of the cryo- 
stat. In terms of dimcnsionless variables, p = kR, fp — 
ep/e+, 7 = TV/R and V = (e+i^)^/^ we have 



dp ^■'p ^ 



(12) 



The solution to this with the boundary condition /|p=i 
1 is 



.n^'=iP-'/'+i-i, (13) 



or, in terms of the dimensional energy spectrum, 



k{kRfl^ 



[kRf/^ 



1 



(14) 



Then, expressing the mean vorticity through this energy 
spectrum, one obtains the closed equation for the param- 
eter 7: 



r = 7i?/F = 9^0(7) 



(15) 



which manifests the existence of several different regimes 
of turbulence. Consider first the case when the resulting 
7 < 1 and Eq. 1)14(1 can be rewritten as 



2„,2 



^2fc5/3 



fc2/3 



where 



v3/2 



(1 -7)3/2i? 



(16) 



(17) 



is a crossover wavenumber separating two different scal- 
ing ranges. For k ^ fccr we have the K41 scaling in 



which the dissipation is negligible and the energy flux is 
approximately constant. 



(18) 



i?2/3fc5/3 ■ 

This equation can be rewritten in the traditional form (jT)) 



Ek 



-00 I 



(19) 



in which the energy flux Eqc for k > fccr due to the mutual 
friction can be much smaller than the energy influx, e+, 
into the turbulent system: 



e+(l -7)3 < £^ 



(20) 



In order to see how Sk decreases toward large fc, ap- 
proaching £00 at fc ~ fccr consider region k <^ fccr- In this 
region Eqs. ifT^ . l(TI|l yield: 



£k = 



{kRf ' 



Ek = 



i?2 fc3 



(21) 



Thus the rate of the energy dissipation, being propor- 
tional to Ek: decreases toward large fc and becomes in- 
significant at fc ^ fccr. In this region Sk — £00 and one 
has the K41 scaling ifT^ . 

The K41 scaling ends by a cutoff at a "microscopic" 
scale 1/fc* at which circulation in the fc*-eddy reaches 
the circulation quantum n: 



RV{l--i) 

(fc*i?)4/3 ' 



I.e. 



(fc,i?)^/3^AA^(l-7), 



(22) 



(23) 



where N"^ = RV/ n is the "quantum Reynolds number" 
(ij . Parameter M can be considered as the ratio of R to 
the mean inter-vortex distance. Clearly, with the clas- 
sical approach to the problem we can consider only the 
limit A/" > 1. 

Now we can clarify the equation of self-consistency H15|l 
for 7. Estimating loq as follows 



= {\^\f (|^|2) 



dk k^ Ek , 



i/R 



and using Eqs. (|15|) and Hlt)|) one gets 



q- I dkk^/^ 

l/R 



fc2/3 



u2/3 



(24) 



Together with Eqs. (|17|l and H23|l . that relate fc* and fccr 
with 7 and TV, this equation allows to find 7, fc, and fccr in 
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the terms of "external parameters" of the problem, q and 
Af. As we pointed out, the classical regime of developed 
turbulence corresponds to the region g <C 1, A/" 1. 

Consider first the case when the inner (quantum) scale 
of turbulence is well separated from the crossover scale: 
A;* S> kcr- Then the main contribution to the integral in 
Eq. (I24II comes from the region k ^ kcr, when the second 
term in the integral dominates. Therefore Eq. (|24(l gives 
the relationship 



(25) 



which together with Eqs. (|17(l and H23|l yields Eq. for 7: 



gA/-^7/(l-7)'/' 



(26) 



For qJV <C 1 the solution is 7 
Rkcr ~ [qAff^ < 1 and: 



Ek 
Rk^ 



k{kRf/^ 



qj\f < 1 that gives 
l~2qJV{l- 

(27) 



(fci?)2/3 

for gA/" < 1 . 



It means that for qj\f ^ 1 in the entire inertial inter- 
val, < fc < fc*, one has usual Kolmogorov-Obukhov 
spectrum with the small, of the order of qj\f, negative cor- 
rections. In other words, at qj\f <^ 1 the mutual friction 
has a negligible effect on the statistics of turbulence. 

The situation is different in the region qM 3> 1 . In this 
case the solution to Eq. is 1 - 7 ~ {qUy^/^ < 1 
and instead of Eq. H27|l one has: 



R^k5/3 



1 



1 



f,2/3 



^2/3 



for 



qj\f, k^R~J\f/q^/^:$>k^,R 



qAf:$>i, 

(28) 



One sees that the pumping and the crossover scales are 
well separated if qAf ^ 1. Under this condition, the 
quantum cutoff scale is also well separated. 

Up to now we assumed that the second term in [. . . ] 
in the integral H24|l is dominant. In the opposite case 
instead of Eqs. H25|) and H2()|l one has: 



k^R : 



kcrR ~ A/"^ exp — - 



(29) 



Taking into account that we are considering solutions 
in which the first term in [. . . ] in the integral (|24(l is 
dominant we have to take fccr > k* which gives Af > 
exp(l/g2). Therefore the range of parameters where the 
two-cascade regime occurs is 



l/q < A/" < exp(l/g2) . 



(30) 



An important feature of this solution is that both the 
energy spectrum Ek and the spectrum of the flow dis- 
sipation Ek are concentrated at the largest length scale 



i?, whereas the dissipation is mediated by the vorticity 
ujo concentrated at the smallest (microscopic) scale 1/fc*. 
Therefore the energy balance between the Kolmogorov 
cascade and the energy dissipation must occur already 
for the largest eddies. This gives the condition 



R 



(31) 



This means that in this turbulent state the mean vorticity 
is 070 = T/q = U/qR. If qN 1 one has 7 close to 1 
and thus our double-cascade solution H28I) satisfies the 
large-scale balance (|31|l . 

At qAf ~ 1 one has fccr — ^/R, i-e. the region of 
the k~^ spectrum shrinks. At qAf 4C 1 the parameter 7 
deviates from unity, 7 ~ qAf. Here two scenarios are pos- 
sible. In the first one we have solution (|27|l in which the 
mutual friction is unessential and thus is unable to com- 
pensate the Kolmogorov cascade. When the intervortex 
distance scale is reached the Kolmogorov energy cascade 
is then transformed to the Kelvin wave cascade Q for 
the isolated vortices. In the second scenario suggested 
in Refs. and [13, at qAf <C 1 the turbulent state is 
completely reconstructed and the so called Vinen state 
emerges. This state introduced by Vinen |1J|] and then 
by Schwarz 0| contains a single scale r = k/V and thus 
no cascade. 

Another interesting case to consider is 7 > 1 when 
Eq. H14f) for Ek can be rewritten as follows: 



Ek 



R^k^/^ 



fc2/3 r2/3 



where 



v3/2 



(7 - l)3/2i? 



(32) 



(33) 



Here the Kolmogorov cascade stops already at the scale 
fccr, and the equation for 7 reads: 



(34) 



The solution satisfying the large-scale energy balance 7 : 
1 is 




7 - 1 ~ exp(-2/3(?2 



(35) 



This solution is self-consistent and does not require the 
microscopic scale cut-ofi' if the circulation at the scale fccr 
is big enough, i.e. Wcr/fccr ^ k. This occurs, however, 
at very large counterfiow TV 3> exp(l/g2). However, we 
think that this solution is unstable. Indeed, let us con- 
sider a perturbation of this spectrum in a form of cutting 
off its tail. This will lead to a reduction in the dissipa- 
tion rate so that the subsequent evolution will build a a 
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Intensity, J\f 


Crossover, fccr^? 


Quantum cutoff, fc*i? 


1 T 

1< AA< - 

1< i < < e^/"' 
1 < e'/"' < 


None, |-scaling 
gAA 

None, 3-scalmg 


AA 



TABLE I: Scaling range boundaries in the cases of weak, in- 
termediate and strong pumping. 



K41 constant-flux tail rather than restore the scal- 
ing. The K41 tail will strengthen until its contribution 
to the friction will restore the energy balance. Thus, the 
resulting new steady state will have the K41 part, i.e. 
will be of the first kind. 

The summary of the different regimes (without possi- 
ble Vinen state) is shown in the Table. 

Now let us compare our results for the superfluid tur- 
bulence with earlier works. In Refs. ^ and the 
effect of the mutual friction on the high momentum tail 
was overestimated, which led to the incorrect result for 
the spectrum of dissipation at high momenta. However, 
some general features suggested in Ref. j5| remain the 
same. In particular, there are different regimes of the 
superfluid turbulence. The transition line between two 
turbulent reg imes, qAf ~ 1, was correctly determined in 
Refs. 0, 0|, as well as the vorticity in the regime of 
cascade, ujq = V/{qR). Recently Vinen used a dif- 
fusion equation model which is similar in spirit to the 
model used in our Letter and, perhaps, even better be- 
cause in properly accounts not only for the cascade states 
but also for the thermodynamic equilibria. However, this 
equation is harder to solve analytically and Vinen used 
numerical simulations, his qualitative and numerical re- 
sults are consistent with our analytical solution. 

In conclusion, we discussed the spectrum of the su- 
perfluid turbulence governed by the nonlinear energy 
dissipation due to the mutual friction between the vor- 
tices and the normal component of the liquid, which re- 
mains at rest. We found that in agreement with Refs. 
P, 1^ , the flow states are determined by two dimension- 
less parameters: the velocity-independent Reynolds num- 
ber Tie = 1/q which separates the laminar and turbulent 
states; and the quantum velocity-dependent parameter 
TV = y^VR/ K which contains the quantum of circulation 
around the quantum vortices k and which determines 
the transition, or crossover, between different regimes 
of superfluid turbulence. In some region of the (q',A/') 
plane, we found the turbulent state with a well defined 
Richardson-type cascade. This state displays both the 
Kolmogorov-Obukhov |-scaling law oc k~^/^ and the 
new "3-scaling law" Ek oc each in a well separated 



range of k. Possible connection of the phase diagram of 
the flow states to experimental observations is discussed 
in Ref. 
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